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Quantum mechanics still provides new unexpected effects when considering 
the transport of energy and information. Models of continuous time quantum 
walks, which implicitly use time-reversal symmetric Hamiltonians, have been 
intensely used to investigate the effectiveness of transport. Here we show how 
breaking time-reversal symmetry in this model can enable directional con- 
trol, enhancement, and suppression of quantum transport. Examples ranging 
from exciton transport to complex networks are presented. This opens new 
prospects for more efficient methods to transport energy and information. 
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Understanding quantum transport is key to developing more robust communication net- 
works, more effective energy transmission, and improved information processing devices. Con- 
tinuous time quantum walks have become a standard model to study and understand quantum 
transport phenomena HI l2l El IS [51 . Time-reversal symmetry (TRS) characterizes all quan- 
tum walk models to date. This symmetry implies that the site-to-site transfer probability at time 
t = T is the same as at time t = —T, thereby prohibiting directional biasing. Here we introduce 
and study continuous time "crural" quantum walks whose dynamics break TRS. Our findings 
show that the breaking of TRS offers the possibility of directional biasing and allows one to 
suppress or enhance transport relative to the standard quantum walk. One reason this effect 
was overlooked could be that time-reversal asymmetry cannot affect the site-to-site transport in 
some simple cases, such as linear chains and trees. This is proven in the supplementary online 
material (SOM) Section 2. Prior efforts have focused on controlling and directing quantum 
transport using either in situ tunable Hamiltonians j6l|71 or tailoring specific initial states (HI. 
In contrast to known approaches, we consider states initially prepared in the standard site-basis 
and time-independent Hamiltonians that induce time-asymmetric evolutions. 

To demonstrate the effect of TRS breaking, we chose four examples which illustrate the 
main ideas of directionality, suppression and enhancement of transport. The first example is 
a unitary quantum switch, where the direction of quantum transport is controlled by a time- 
reversal asymmetry parameter. The second example examines transport in a linear chain of 
triangles, showing a 633% transport speedup for the chiral quantum walk. Next, we demon- 
strate complete suppression of chiral quantum walks on loops with an even number of sites. We 
then consider a system widely studied in the exciton transport literature: the Fenna-Matthew- 
Olsen complex (FMO). Although this naturally occurring system is highly efficient, we find that 
the introduction of chiral terms allows for an enhancement of transport speed by 7.68%. Finally, 
to investigate the robustness of TRS breaking effects on transport, we consider randomly gen- 
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erated small-world networks. By appending time-reversal asymmetric terms to only the edges 
of the graph connected to the final site, we could increase the speed of the site-to-site transport 
significantly, up to 130%. 

In the standard literature on continuous time quantum walks lHJ|2]|3l|4l[5]|, the time-independent 
walk Hamiltonian is defined by a real weighted adjacency matrix J of an underlying undirected 
graph, 

sites 

H QW = ^ J nm (\n) (m\ + \m) (n\) . (1) 

The condition that the hopping weights J nm are real numbers implies that the induced transitions 
between two sites are symmetric under time inversion. We can break this symmetry while 
maintaining the hermicity of the operator by appending a complex phase to an edge: J nm — > 
Jnmfi %6nm resulting in a continuous time chiral quantum walk (CQW) governed by 

h cqw = J nme Wnm \n) (m\ + J nm e~ ienm \m)(n\. (2) 

n,m 

When acting on the single exciton subspace the Hamiltonian given in Eq. © can be expressed 
in terms of the spin-half Pauli matrices: 

H C QW = ^ Jnm COS ( 6 nm) + ^m) + Jnm Sin ( 6 nm) - 0-n°m) ( 3 ) 

n,m n,m 

which arises in a variety of physical systems. In the CQW framework, we investigate coherent 
quantum dynamics and incoherent dynamics within the Markov approximation. Both types of 
evolution are included in the Lindblad equation lf9l ITOl ITTL IT2]] : 

—p(t) = £{p} = -i[H CQ w, p] + Y2 Lk P L k ~ \ ( L k L kP + pL\L k ) (4) 

k 

where p(t) is the density operator describing the state of the system at time t and L k are Lindblad 
operators inducing stochastic jumps between quantum states. For example, using the usual 
terminology of Markovian processes, we call site t a trap if it is coupled to site s by the Lindblad 
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jump operators, L k = \t)(s\. The site-to-site transfer probability, P n _ vm (t) = (m\p(t)\m), gives 
the occupancy probability of site m at time t with initial condition p(0) = |n)(n|. To quantify 
the transport properties of quantum walks, we use the half-arrival time, n/2, as the earliest time 
when the occupancy probability of the target site is one half. We will also make use of the 
transport speed, uy 2 , defined as the reciprocal of n/2- 

We now introduce a quantum switch which enables directed transport and could, in princi- 
ple, be used to create a logic gate and offer future implementations of transport devices to store 
and process energy and information. Fig. Q] presents an example of this switch. The value of 
a phase (e td ) appended to a single control edge across the junction allows selective biasing of 
transport through the switch. The maximal biasing occurs at \9\ = n/2, and the sign determines 
the direction. The first maxima of P s ^ E (t) (transfer probability from site S to E) in the unitary 
dynamics without traps can be enhanced by 134% or suppressed to 91% with respect to the 
non-chiral case. When considering traps in the Lindbladian evolution, the optimal transport ef- 
ficiency is 81.4% in the preferred direction. The switch violates TRS as P s ^e(— t) ^ Ps~^e{P)- 
By using P s ^E(—t) = PE-*s(t) and the symmetry of the configuration P E ~>s(t) = Ps-*F(t), 
we conclude that transport is biased towards the opposite pole when running backwards in time, 
see Fig [TJ Note that the behavior of the switch is largely independent of the length of the con- 
necting wires. 

We will now utilize the directional biasing of the triangle to give an example of a speedup 
of chiral walks. Using the composition of eight triangular switches as depicted in Fig. I2A[ by 
simultaneously varying all phases along the red control edges to the same value, we examine 
the effect of time -reversal asymmetry on transport. We find that the occupation probability as 
a function of 9 is symmetric about ±tt/2 with the negative value corresponding to maximal 
enhancement and the positive value to maximal suppression. Unlike the occupation probability 
maxima in the switch, here the first apexes are separated in time. When we include trapping, 
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Figure 1: The quantum switch. (A) Directional biasing: enhanced transport in the preferred 
direction. The plot shows the occupancy probability Ps^e of site E with the particle initially 
starting from site S with and without sink (dashed and solid lines, respectively). (C) This 
evolution is time-reversal asymmetric as replacing t with —t results in the particle moving from 
site S towards site F. (D) When starting at site E, the particle evolves towards site F. (E) By 
replacing t with —t a particle initially at site E evolves towards the initial configuration (B). To 
recover time-reversal symmetric transition probabilities in the evolution (B), requires that one 
also performs the antiunitary operation lfT3ll on the Hamiltonian mapping 9 to —9. This has the 
same effect as reflecting the configuration (C) horizontally across the page while leaving the 
site labels intact. 
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Figure 2: (A) Triangle chain and (B) the FMO complex. (A) The phase e is applied to the red 
edges simultaneously in the triangle chain. The plot illustrates the occupancy probability at the 
end site E as a function of time for different values of the phase 6 with and without trapping 
(dashed and solid lines, respectively). (B) shows the occupancy difference with respect to the 
time reversal symmetric Hamiltonian of the FMO complex. We use an optimization procedure 
to enhance the transport. While holding the magnitude of the couplings constant, we optimize 
two sets of phases, A\ and A 2 , which correspond to seven and three edges with an enhancement 
at 7~i/2 of 3.25% and 2.25%, respectively. 
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the half-arrival time is reduced from the non-chiral value Ti/ 2 = 38.1 to 5.2 which represents 
a 633% enhancement. To conclude this section we focus on suppression of transport by chiral 
quantum walks. A good example is the polygon with an even number of sites. In this case 
complete suppression can be achieved by appending a phase of ir to one of the links in the cycle; 
thereby rendering it impossible for the quantum walker to move to the diametrically opposite 
site. The proof that this site-to-site transfer probability is zero for all time can be found in SOM 
SI. 

In natural and synthetic excitonic networks such as photo synthetic complexes and solar 
cells, we are faced with non-unitary quantum evolution due to dissipative and decoherent inter- 
action with the environment. Studies have shown that dissipative quantum evolution surpasses 
both classical and purely quantum transport (for interesting recent examples see [IHEHl). A 
widely studied process of such dissipative exciton transport is the one occurring in the Fenna- 
Matthews-Olsen complex (FMO), which connects the photosynthetic antenna to a reaction cen- 
ter in green sulfur bacteria 021 QUI. Due to the low light exposure of these bacteria, there is 
evolutionary pressure to optimize exciton transport. Therefore, the site energies and site-to-site 
couplings in the system are optimized, yielding a highly efficient transport. However, the natu- 
ral exciton transport occurs via time-reversal symmetric interactions, and we ask whether TRS 
breaking interactions may further enhance the efficiency of the process. We consider the real 
Hamiltonian modeling transport on the FMO and allow for TRS breaking by introducing com- 
plex phases and find that the transport speed can be further increased. We study the seven site 
model of the FMO using an open system description that includes the thermal bath, trapping at 
the reaction center, and recombination of the exciton [15 J. By performing a standard optimiza- 
tion procedure (as outlined in the SOM S3) that varies the phase on a subset of seven edges, 
we found a combination of phases where the transport speed, u\/2, is enhanced by 7.68%. In 
Fig. I2BI the enhancement of the time dependent occupation probability is shown for the chiral 
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quantum walk. We note that optimization over only three edges already changes the transport 
speed by 5.92%, see SOM S3. 

One may also ask if optimization procedures limited to small subsets of edges will generally 
lead to improved transport in larger and possibly randomly generated networks. We found a 
positive answer when testing the site-to-site transport between oppositely aligned nodes in the 
Watts-Strogatz model ifTTl . This family of small- world networks continuously connects a class 
of regular cyclic graphs to that of completely random networks (Erdos-Renyi models lfT8lO by 
changing the value of the rewiring probability. 

We numerically investigated graphs with 32 nodes, average degree 4 and range over rewiring 
probability p considering 200 different graph realizations for each value of p. An example 
with p = 0.2 is depicted in Fig. 13 Al Here the occupancy of a sink connected to site E is 
compared between the chiral walk and its achiral counterpart. The particle begins at site S and 
we perform the optimization of the phases only on edges connected to site E. In the case of the 
chiral quantum walk, the sink reaches half-occupancy in 54.8% less time on average. When the 
occupancy of the sink of the standard quantum walk reaches 0.5 at time T = 175.17, the sink 
occupancy of the chiral walk is already 33.5% higher on average. 

In all the examples studied, we found that transport is affected by the topology of the net- 
work. A key explanation is the following. If two Hamiltonians are related by on-site unitary 
transformations mapping |n) to e ia " \n), then the site-to-site transition probabilities will be iden- 
tical. This fact provides a tool to reduce the effective space of phase parameters for controlling 
transport, see SOM SI. For instance, the site-to-site transfer probability is insensitive to phases 
in tree graphs. For bipartite graphs the phases can have an effect, however, the dynamics still 
remains time-reversal symmetric. 

A further consequence is that the sums of phases along a chosen orientation of a loop are 
the unique invariants under on-site unitary transformations. In a wide range of cases and par- 
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Figure 3: Transport enhancement of the chiral quantum walk is robust across randomly gener- 
ated Watts-Strogatz networks. An example of this small-world network, with rewiring proba- 
bility p = 0.2, is depicted in (A). The transfer probability P from site S to the sink connected 
to site E is plotted in a realization of the network. (B) shows the average enhancement of half 
arrival time (Ati/ 2 ) for different values of p. 
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ticularly in all examples we considered, we found strong evidence of the robustness the effect 
has on transport. For instance, the examples in Fig. I3BI shows that in the Watts-Strogatz model, 
the transport enhancement due to the time reversal asymmetry of the Hamiltonian is insensitive 
to changes of the rewiring probability p and the clustering coefficient measuring the density of 
triangles in the graph. Finally, additional calculations show that scale free networks such as the 
Barabasi-Albert model [fT9l , show a similar transport enhancement, indicating robustness also 
with respect to the degree distribution. 

This study pioneers the exploration of a new degree of freedom that allows for a significant 
improvement of control in the engineering of quantum transport. The fact that we were able 
to optimize and control transport by adjusting the phase on only a few edges inside a complex 
network and that the effect was relevant in a host of examples adds optimism to the robustness 
of this approach. Experimental demonstrations of the effects we predict are within reach of 
existing hardware, as outlined in SOM Section S2. 
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Supplementary Material 

SI Basic definitions and properties of chiral quantum walks 

51.1 Site-to-site transfer probability 

The Markovian open-system dynamics of a continuous time chiral quantum walk is given by 
the Kossakowski-Lindblad equation 191 fIDl [TTl [121 

£{P} = -i[HcQW,P\ + ^2 ° mn \ L ™P L nm ~ ^nm^nmP ~ ljP^\ m ^nra \ , (SI) 

(n,m) 

where the chiral Hamiltonian Hqqw is defined in Eq. (2), and the Lindblad operators are given 
as L mn = \m) (n\ with c nm > 0. Transport from vertex IS") to vertex \E) during such dynamics 
is characterized by the site-to-site transfer probability (STP). In the unitary case (c nm = 0) it is 
given by 

Ps^ E (t) = Tr (e-^^ps e iHc ^ p E ) (S2) 
with p s = \S)(S\ and p E = \E)(E\, while for the general Markovian case it is 

Ps^ E (t) = Tv(e ct {p s }p E ) . (S3) 

51.2 Time-reversal symmetry of the unitary achiral dynamics 

In the setting of quantum walks, the time-reversal operator T acts as complex conjugation (with 
respect to the vertex basis) lfT3l : 

T^a» = ^2a* v \v) . 

The antiunitarity of T and T 2 = 1 implies that = T. The time-reversal of a Hamiltonian H 
is given as THT\= THT). The H >-)■ THT action is represented in parameter space by the 
replacement 9 mn i->- — 9 mn in Eq. (2). Thus exactly the achiral quantum walks are left invariant 
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by this action. The STP's of H (P s ^ E (t)) and that of H' = THT (P' s ^ E (t)) are related in the 
following way: 

P's^ E (t) = Ps^E(-t) and P' s ^ E (t) = P E ^ s (t), 
which can be verified using Tp v T = p v and the cyclicity of the trace as follows: 

Ps^ E (t) = Me^ THT »pse^ HT »p E ) = Tr(Te im Tp s Te- im Tp E ) 

= Tr(e iHt Tp s Te- iHt Tp E T) = Tr(e lHt p s e~ im p E ) = P s ^ E (-t) , 
Ps^E(-t) = Tx{e iHt Ps e- mt p E ) = Tr(e~ lHt p E e lHt Ps ) = P E ^s(t) • 

A crucial consequence of the above is that in the case of achiral quantum walks, the transition 
probabilities are the same at time t and —t, i.e. Ps-+ E (t) = Ps-> E (—t), and directional biasing 
is prohibited P s _> E (t) = P E _> s (t). However, H ^ THT^ does not necessarily imply that 
transition rates are asymmetric in time. This is because THT^ might be gauge-equivalent to H, 
as will be seen in the next section. 

SI. 3 Gauge transformations 

A local change of basis, such as a diagonal unitary: 

U d \n) = e ian \n) , (S4) 

leaves the STP invariant. To prove this, let us first note that any unitary basis-change U would 
induce a transformation on the Lindblad superoperator £ — > £ with 

£'{p} = U£{U ] P U}U^ . 

Using this and the invariance of localized states under diagonal unitaries (U\p v U d = p v ), we 
arrive at 

P's^eH) = Tr(e c,t { Ps } PE ) = Tv{U d e a {Ulp s U d }Ulp E ) 

= Tx{e a {p s }Ulp E U d ) = Tr(e a {p s }p E ) = P s ^ E (t) , 
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which proves the invariance of the STP under the gauge transformations defined by Eq. (|S4T) . 

Under these diagonal transformations, the parameters of the quantum walk Hamiltonian 
transform as 

@mn ' ^ @mn &n j (S5) 

as illustrated in Fig. |ST] The incoherent part of the Lindblad equation dSTI) does not change 
since the Lindblad operators transform as L nm — > e^ am ~ an ^ L nm and these phases cancel in 
Eq. dSTI) . since L nm and L\ m appear paired. Two important properties of the model now follow: 
(%) phases on tree graphs can be transformed out completely and (ii) the sum of phases along 
loops is invariant under gauge transformations. 




Fig. SI: Effect of the gauge transformation \n) — > e ia " \n) on vertex n. Phases on edges can be 
gauge-transformed without changing the transition amplitudes, as described in the text. Here 
we arrange the graph as a tree rooted at d — 0. 

The first property is illustrated in Fig. [STJ Let us take an arbitrary tree graph and pick a 
vertex m with only one neighbor. Redraw every other vertex on successive levels characterized 
by the distance d of the vertexes from the given vertex m. Note that the number of edges 
connecting two vertices, d, is by definition, unique in tree graphs. In such an arrangement only 
one edge emanates downwards from a given vertex on a line of d > so Fig. ISTl represents 
the general neighborhood of a vertex n having distance d = 1 from m. The indicated gauge 
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transformation with a n = —9 removes the phase from the bottom edge. Then, one iterates 
the procedure for all vertices at level d = 2 and consecutively for all levels. In this way, all 
phases are removed. For the second property, pick an orientation on a loop of N vertexes and 
compute 9 := J2iLi considering 4>n,n+i = </>n,i- A gauge transformation |n) !->■ e ia "\n), 

according to Eq. (IS5l) leads to: 



so the sum remains unaffected. 



SI. 4 Analytic examples for enhancement, suppression and direction 

Using the example of the regular polygon (that has J mn = 1 for all edges) we will explain 
how enhancement, suppression and direction of transport can be influenced by the phases in 
chiral quantum walks. Typically we need to solve the eigenvalue problem of the Hamiltonian 
H\j) = Ej\j) (the tilde denotes energy eigenstates) to compute 



S-^E 



(t) 



N 



J>lM|S->e- 



For a regular homogeneous polygon (H = ^2^=1 & lip \ n ){ n + 1| + h.c. with cyclic boundary 
condition) the eigenvectors and eigenvalues are found from the Fourier Transform 



1 N 



e N 

The STP from site S to E then reads 
We now list several features 



\n) , E k (ip) = 2 cos 



2nk 



N-l 



1 \- -i^cce^-vH^^ 

N f-^ 



k=0 



(S6) 



Direction and enhancement 
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In the case of time symmetric quantum walks, for the homogeneous polygon Ps->e = 
Ps^n-e+2S- This implies that if the number of edges is odd, the STP is at most 1/2, since 
there is always a vertex different from E at which the probability is identical. The chiral 
case breaks this symmetry. In particular, for the regular triangle one can reach P 1 _^. 2 ('t) = 
1; one can find values of 9 analytically by equating the phases of the three exponentials 
in the STP modulo 2n and solve the two linear equations for t and 9. Enhancement and 
direction is indicated in Fig. |S2l where a triangle graph is shown with inhomogeneous 
coupling. 

• Complete suppression of transport 

Complete suppression of transport is possible in loops with an even number of sites. They 
belong to the set of bipartite graphs, whose vertices can be partitioned into two sets such 
that edges connect only members from different subsets. Directing transport in bipartite 
graphs is forbidden. This can be easily shown by noting that the gauge transformation 
\n) i — v —\n) for all members of one subset has the effect of H i-> —H, which is precisely 
the time reversal transformation. Hence 

P s ^ E (t) = Ps^E(-t) = Tr(e* Ht p s e-* Ht p E ) = Tr(e- iHt p E e lHt p s ) = P E ^s(t) . 

Even chiral loops suppress transport as 9 = n results in Pi^n/2+i (t) = for all times 
t. This can again be shown by writing down the formula (|S6|) and separating it into 
two sums, one for k being even and the other for k odd. The simplest even polygon 
is the quadrilateral that can be used as a building block for the realization of transport 
suppressing topologies. 
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Fig. S2: (A) Triangle with inhomogeneous coupling and chiral phases. (B) Illustrates the quan- 
tity Pi_t 2 (t) for the triangle graph with inhomogeneous coupling J 12 = 1, J23 = 1.3, J13 = 0.5. 
Directing the phases of the chiral walker effects the transport, in particular, there is enhance- 
ment for 9 = 3ip = n / 2 in both arrival time and maximal probability. (C) The probability varies 
as a function of the chiral angle at a given time near the first maximum. 
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Fig. S3: Half arrival time scaling. Dependence of the half arrival time, n/2, on the length of the 
triangle chain shown in Figure [2K for different values of the phases. Both chiral (9 ^ 0) and 
achiral (9 = 0) cases exhibit linear scaling behaviour. The slope can be tuned by changing the 
value of the phase parameter. 

SI. 5 Robustness of chiral transport enhancement of the triangle chain 



In the triangle chain discussed in the main text, the half arrival time n/2 depends linearly on the 
size of the system, see Figure [S3J This fact is a characteristic of quantum walks in linear chains 
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and continues to hold for the chiral triangle chain. 

S2 Experimental proposal 

52.1 Summary 

In this section we outline a proposal to simulate chiral quantum walks in a system of ultra- 
cold trapped atomic ions. The proposal can be realized with currently available techniques and 
technology. Trapped ion systems have previously been used to investigate non-chiral quantum 
walks [|20l . where the motional state was used to encode the walk. In contrast, in our proposal 
both chiral and non-chiral walks are encoded into the electronic state of the ions. The walk 
dynamics is generated by ion-ion interactions mediated by joint vibrational modes of the ion 
string. These interactions can be driven by laser-induced optical dipole forces, for example, 
and can be implemented with high quality as shown in several recent works ||2T1 122| . In the 
text below we first introduce the 3-ion system that the experiment requires. Then we explain 
how to engineer the quantum walk Hamiltonians that we wish to investigate. Next, we calculate 
the chiral and non-chiral phenomena that can be generated by these Hamiltonians. Finally, we 
outline the experimental procedure and note some important experimental considerations. 

52.2 Implementation in trapped ions 

We consider a string of three ions in a linear ion crystal, which can be achieved with standard 
linear- Paul traps. A long-lived electronic transition inside each ion encodes a two-level spin. 
The mapping could be that spin-down is represented by an ion in the ground state, while spin- 
up is represented by an ion in the excited state of this transition, for example. The 8 possible 

logical states of this system then become | |, I, !), | !, t)> I > I T>T>T}- In our 

proposal a subset of these states represents the sites of the quantum walk, as will be described. 
The dynamics of any quantum walk between these states require the ability to turn on interac- 
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tions between the ionic spins. Effective ion-ion interactions can be achieved by off-resonantly 
coupling the electronic states in which the spin is encoded to one or more common vibrational 
modes of the string. Such a spin-dependant force can be implemented using a laser field with 
two symmetrically tuned frequencies at cu s ± \i with fi <C u s . If the symmetric detuning is suf- 
ficiently far from all motional sidebands so that the generation of phonons can be adiabatically 
eliminated [|23l . then a pure spin-spin interaction is generated of the form: 



Here describes the laser coupling strength to the electronic transition, r\i is the Lamb-Dicke 
parameter describing the coupling strength between the laser and the motional sideband for a 
single ion. An experimentally tunable phase 4> is set by the phase difference of the two driving 
laser fields. Finally, describes how strongly ion i couples to motional mode u. The increasing 
difficulty of generating spin-spin interactions with more ions in a string is incorporated in 

In this proposal we will simultaneously use two of the normal modes along the axis of 
the string, called the center-of-mass (COM) and breathing (Br) modes, to mediate spin-spin 
interactions. These modes are typically well- spaced in frequency and can therefore can be 
individually addressed. The coupling vectors fej for the axial COM and breathing are 1 / \/3 x 
[1, 1, 1] and l/\/2 x [1, 0, —1], respectively. Consequently, the contributions to the spin-spin 
interaction Hamiltonian of these modes are: 




(S7) 



i<3 



where 




y 



(S8) 



Hcom{4>i) 



t I '1 I 2 3 | 1 3 \ 

JCOM (0-0! 0-0! + 0^0^ + ^!^J 



(S9) 



H Br ((h) 



J Br CTfoVfo 



(S10) 
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where 

A — bj 



determines the coupling strength of the spin-spin interaction due to vibrational mode i. This 
equation allows for separate pairs of light fields to drive each mode simultaneously, which can 
therefore have different coupling strengths and detunings. This is required because we wish 
to be able to precisely control the relative strengths Ai, which can be achieve by changing the 
relative strengths of the laser fields. Separate pairs of light fields are also necessary since we 
require the ability to individually set the phase fa determined by the phase difference between 
the laser fields in each pair. 

S2.3 Ideal experimental behavior of the proposed walks 

Now that we have established the Hamiltonians to drive the quantum walks, given by equations 
(|S9l) and (IS 101) , we will examine the walk dynamics that we wish to explore. A key observation 
is that the operators JaTa^ (m, n £ {1, 2, 3}, m ^ n) leaves invariant the subspace spanned 
by the following four spin states 

|1> = | t,l,l>, |2> = | |,t,l), |3> = | U,t), |4> = | t,t,t> • 

In order to see this, consider the following decomposition (m, n £ {1,2, 3}, m ^ n): 

a™al = (cos 0a™ + sin 0a™) (cos 0< + sin 0^) (SI 2) 

in terms of the creation and annihilation operators a± = a™ ± ia™. It follows that: 

a <f> a <f> = e ° + a + ~| e °- a - + (T + U_ + cr_u + (olj) 

Here \n) and \m) are coupled by the real strength J, whereas |4) is coupled to the remaining site 
by the same strength multiplied with the phase e~ 2t ^ . Now, notice that the modes H C om and 
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Fig. S4: Networks describing the modes H Br (left) and H C om (right). 

Hbt are built from operators of the above form. Using these models we can realize quantum 
walks on the four sites given above. Now one determines the parameters of the quantum walk 
Hamiltonians: 

4 

# (fc) = E E JffiWH with Jffi=^. 

k n,m=l 

where the index k refers to the two different modes Hcom{4>i) and H B r{4>2)- Figure IS4l shows 
the general situation. 

The values to be used for the experiment are given in the following table: 





J COM 


J Br 


01 


02 


HcQWx 


2 


-3 


vr/2 


0.304 n 


Hcqw 2 


2 


-3 


tt/2 


-0.304 7T 


Hqw 


-2 


1 


tt/2 






The parameter 2 is determined by cos2</> 2 = —1/3. Since the effect 0i is just a minus 
sign in Room for the edges incident to |4), we see that Hcqwx and H C qw 2 are the images 
of each other under TRS transformation. One can also check that i?Qw/ is the time-reversal 
symmetric counterpart of the chiral walks in that all its couplings are real and equal to the 
absolute values of the chiral walks after we relabel sites |2) and |4). Computer simulation 
(Figure [S5]) shows the effect of TRS breaking, since P 1 _ >2 (*) P2->i(t), enhancement and 
suppression of transport as the first maxima of Hcqw 2 an d Hcqw! are of considerably larger 
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Fig. S5: The effect of the chiral phase parameter on the occupation probability. The plot shows 
the chiral walkers at site 2 and the achiral one at site 4, all being initially at site 1. By this 
arrangement the real coupling strengths of the corresponding Hamiltonians are the same. 

(smaller) magnitude, respectively, than that of the corresponding symmetric walker. When we 
compute P 2 ^i(t), the roles of the two chiral walks are exchanged due to the symmetry of the 
model. 

S2.4 Experimental details 

As a specific example we consider using three 40 Ca + ions in a standard linear Paul trap. This 
atomic species is used by several groups around the world and can be precisely manipulated in 
the way that we require. Spins can be encoded into a metastable electric quadrupole transition, 
the excited state of which has a lifetime of approximately 1 second which is much longer than 
the times required for quantum operations. Spin operations can be driven using a 729nm laser, 
spin-spin interactions can be implemented using a bichromatic light field with symmetrical 
detuned sidebands around the electronic transition 11241 . We note that implementation of our 
proposal using another atomic species with a hyperfine transition, for example, would also be 
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possible. 

FigurelS6lshows one possibility allowing simultaneous realisation of the Hamiltonians given 
in equations (|S9l ) and (IS 101) . Two pairs of bichromatic fields are turned on, on pair tuned close 
to the axial COM and another to the axial breathing mode. We choose to use the axial modes, 
since they are well- spaced out in frequency, thereby enabling each bichromatic light field to be 
simultaneously close enough to its nearest vibrational mode such that the effect of that mode 
completely dominates the spin- spin interaction and sufficiently far away such that the adiabatic 
approximation still applies. 

The criterion to generate a pure spin-spin interaction with each bichromatic field is \co m — 
lA 3> Vm^m, the detuning from all sidebands is much greater than the coupling strength 
on those modes. Using the detunings shown in Figure [S6l and {1=2tt x IOOKHz we obtain 
ratios \oj m — fi\/rj m u m of 21 and 14 for the COM and breathing bichromats, respectively. In 
this regime, errors due to the generation of phonons are extremely small. Regarding the issue 
of off-resonant coupling to unwanted modes: straightforward calculations using equation (IS lib 
and the frequencies shown in figure |S6] show that in both cases the size of the far-offresonant 
spin-spin coupling, due to bichromatic field 2 on the COM for example, is more than one order 
of magnitude less than the desired couplings in every case. This could be further reduced by 
increasing the frequency separation between axial modes by increasing the confining potential 
in this direction, for example. 

An important experimental consideration is the requirement to maintain a fixed phase re- 
lationship between the two pairs of light fields generating the walk Hamiltonians. For exam- 
ple, this can be achieved by generating them in the same acousto-optic modulator (AOM) and 
thereby keeping their paths common mode, between the point of generation and the point of 
interaction with the ions. The maximum frequency splitting is approximately 3.4 MHz, which 
will result in some angular divergence of the different frequencies at the output facet of the 
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Fig. S6: Mode frequency spectrum of a linear three ion crystal of 40 Ca + in a typical har- 
monic trap with radial COM (not shown) and axial COM frequencies of 5 MHz and 1MHz, 
respectively. Two bichromatic laser fields, each symmetrically detuned from the axial COM 
(by lOOKHz) and breathing modes (by 50KHz), drive interactions between spins encoded into 
an electronic transition of each ion. The phase difference between the spectral components in 
bichromatic field 1 (2) determines the phase 4> in the Hamiltonian given in equation (|S9l ) ( (1S10I) ). 

AOM crystal, which can be compensated with linear optics allowing coupling into an optical 
fibre for transport to the ion trap itself. 

The experiment would proceed as follows. Firstly, standard methods of doppler cooling, re- 
solved sideband cooling on the axial COM and breathing modes, and optical pumping prepares 
the three ion string into an ultra low entropy state and the initial spin state | I, I, \). Next, a 
standard combination of single-ion focused and three-ion focused beams is used to prepare the 
initial state | f, |, 1). The bichromatic light fields which implement the desired walk dynam- 
ics are then turned on for a fixed period of time. Finally the state of each encoded spin (up 
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or down) is measured individually using standard fluorescent detection techniques and a CCD 
camera. Experiments are repeated multiple times from which estimates of the probability for 
finding spins in all configurations can be deduced. 

S3 Fenna-Matthews-Olson complex 

The Fenna-Matthews-Olson (FMO) complex is a well studied light harvesting system in green 
sulfur bacteria. The interest in FMO resides in the high efficiency shown in the transport of 
energy from the antenna to the reaction center. The faithful description and simulation of the 
FMO complex has been the highlight of recent research. 

The system is composed of three subunits, each made of seven BChl-a molecules embed- 
ded in a protein scaffold. We restrict our study to the simulation of the exciton transport in a 
single unit. To simulate the evolution of the system in the one-excitation manifold, we use the 
following Liouville equation: 

p = C[p] = -i[H, p] + fL^p] + 1 L J {p} + rL T {p} (S14) 
where the Lindblad super-operators are defined as: 

L<f>, n = \n)(n\ (S15) 

L 1>n =\d){n\ (S16) 

L T , n = \r){n\ (S17) 

where the sums are over the site basis and T is the set of sites connected to the reaction center 
(in this case only site three). The coherent part of the evolution is described using a time reversal 
symmetric Hamiltonian from the literature E51 . In (IS 141) the non-coherent terms describe the 
coupling of the system with a thermal bath, dephasing, Eq. (IS 15b , and recombination, Eq. (IS 1 6b . 
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L 4>{P) = Yl L <t>^P L \,n ~ j ( L l,n L <l>,nP + P L \,n L 4>^ 
n 

L y {p} =Y L i,npL\^ n - - (L\ n L j7n p + pL\ n L 1)n ) 

n 

L Ap} = L T) „pLj. )n - - {L\ n L T , n p + pL\ n L T , n ) 
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variance 


A 2 


variance 


3-4 


1.31484899 tt 


4.33 • 10 6 tt 


1.58371001 tt 


4.86 • 10 6 tt 


4-5 


1.66997830 tt 


5.20 • 10 6 vr 


1.39551582 tt 


6.60 ■ 10 6 vr 


6-7 


1.8406103 tt 


1.25 • 10 5 vr 


0.1338368 tt 


1.72 ■ 10 5 vr 


2-7 


1.2949616 tt 


2.80 • 10 5 vr 






1-6 


1.67543320 tt 


5.36 • 10 6 vr 






1-3 


0.04222214 tt 


4.36 • 10 6 tt 






3-6 


0.8761298 tt 


1.03 • 10 5 7T 







Table SI: Transport enhancement in the FMO complex. Both for A\ and A 2 , complex phases 
are applied to the red edges. The results of a simultaneous optimization procedure are listed in 
the table. 

and the effect of the excitation trapping at the site connected to the reaction center, Eq. (IS 17b . 
The thermal bath, for the present paper, is considered as a set of harmonic oscillators coupled 
to the system of interest as in [|T5l ITTTl . 

The initial state is set on site one and the reaction center is connected as an energy sink to 
site three. In the simulation we use a dephasing rate of <\> = 9.0 ps" 1 (which corresponds to 
295K), the recombination rate is 7 = 1 ns~ x and trapping rate from site three to the reaction 
center is r = 1.0 ps -1 . 

The optimization of the phases was found to be robust with respect to phase changes. The 
resulting phases are reported in Table ISTl 
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S4 Small-world networks 



The Watts-Strogatz model ifTTl gives a constructive algorithm for building a network with small- 
world characteristics, starting from a regular lattice. The latter is defined as a periodic chain of 
N nodes where each node is connected to k neighbors (k/2 to the left and k/2 to the right). The 
final small- world network is obtained by taking all neighboring edges at each node and rewiring 
all the edges toward the node left with a probability p. The limit for p — > 1 of the Watts-Strogatz 
model is the Erdos-Renyi random model with fixed number of edges (restricted to connected 
graphs). 

In our case, we set N = 32 and k = 4 and choose only connected graphs as we are 
concerned in comparing transport on networks of the same sites. The initial state is on site S 
while the sink is a external site connected to site E. Site E latter is placed on the opposite 
side of the initial circle. The absorption rate of the sink is r = 1.0. The evolution is described 
by the coherent Hamiltonian, which, in this case, corresponds to the connectivity matrix of the 
network. We define t®Y me ti me needed by the normal quantum walker to reach the probability 
Vsinkif) = 1/2- We add phases to the edges neighboring node E and optimize them in order 
to improve the probability of being trapped in the target site at time Tw 2 . The non-coherent 
part (trapping) is described by a Lindbladian super-operator. We repeat the evolution for 200 
different realizations of the graph for each value of p. 

S5 Numerical Methods 

We used the Quantum Information Toolkit, primarily developed by Ville Bergholm. This is 
a toolkit for the Matlab programming language. The optimization procedure used in the FMO 
and the Watts-Strogatz examples rely on the Interior-point Optimization algorithm of the Matlab 
minimization toolkit. We start the optimization procedure several times from different randomly 
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chosen points of the parameter space, to reach the global minimum of cost function with more 
certainty. Source code for all simulations done in this work is available upon request. 
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